SMA observations of Infrared Dark Clouds: A tale of two cores 
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ABSTRACT 

We present high-angular resolution sub-millimeter continuum images 
and molecular line spectra obtained with the Submillimeter Array to- 
ward two massive cores that lie within Infrared Dark Clouds; one ac- 
tively star-forming (G034. 43-1-00. 24 MMl) and the other more quiescent 
(G028. 53— 00.25 MMl). The high-angular resolution sub-millimeter continuum 
image of G034.43-f00.24 MMl reveals a compact (~ 0.03 pc) and massive (~ 
29 Mq) structure while the molecular line spectrum shows emission from nu- 
merous complex molecules. Such a rich molecular line spectrum from a compact 
region indicates that G034. 43-1-00. 24 MMl contains a hot molecular core, an early 
stage in the formation of a high-mass protostar. Moreover, the velocity structure 
of its ^^CO (3-2) emission indicates that this BO protostar may be surrounded 
by a rotating circumstellar envelope. In contrast, the sub-millimeter continuum 
image of G028. 53— 00.25 MMl reveals three compact (< 0.06 pc), massive (9- 
21 Mq) condensations but with no lines detected in its spectrum. We suggest 
that the core G028. 53— 00.25 MMl is in a very early stage in the high-mass star- 
formation process; its size and mass are sufficient to form at least one high-mass 
star, yet it shows no signs of localized heating. Because the combination of high 
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velocity line wings with a large IR-mm bolometric luminosity (~ 10^ L©) indi- 
cates that this core has already begun to form accreting protostars, we speculate 
that the condensations may be in the early phase of accretion and may eventually 
become high-mass protostars. We, therefore, have found the possible existence 
of two high-mass star-forming cores; one in a very early phase of star-formation 
and one in the later hot core phase. Together the properties of these two cores 
support the idea that the earliest stages of high-mass star-formation occur within 
IRDCs. 

Subject headings: dust, extinction-stars:formation-ISM:clouds-infrared:ISM-radio 
hnes:ISM 



Introduction 



Emerging evidence suggests that the very earliest stages of high-mass star and cluster 
formation occur w i thin cold, dense molecular c lumps called infrared dark clouds (IRDCs; 



Simon et al.ll2006bl : iRathborne et al.ll2006l . 120071 ). Because these molecular clumps have very 



high column densities 



1023-1025 cm-21 



and low temperatures (< 25 K), they have predom- 



Egan et al. 


1998; 


Carev et al. 


1998, 


2000l: 


Simon et al. 


2006a) 



are now un- 



dergoing the process of fragmentation and condensation, IRDCs are important laboratories 
to study the pristine, undisturbed physical conditions of cluster-forming clouds before they 
are shredded apart by stellar winds and radiation. 

Millimeter /sub-mil l imeter continuum studies of IRDCs ( Lis fc Carlstromlll994 Carev et al, 



2000l : ICarav et al.ll2004J : lOrmel et al.ll2005l : iBeuther et al.ll2005l : IRathborne et al. 



2005 



200 



reveal that IRDCs contain many compact cores. These cores have typical sizes of < 0.5 pc 
and masses of ~ 120 ( iRathborne et al.ll2006l ). While most of these cores contain little 
evidence for active star formation, some are associated with bright 24 /xm emission, broad 
molecular line emi ssion, shocked gas, and maser emission, indicating that they ar e actively 



forming stars (e.g.. IRathborne et al.ll2005l : IWang et aL 



a number of low- and intermediate-r aass protostars (ICarev et al 



in addition to high-mas s protostars (jBeuther et al.ll2005l : IRathborne et al 



2006 



Chambers et al. 2008). Indeed 



2000 



Redman et al 



2005 



200 



Pillai et al. 



20061 : IWang et al.ll2006l ) have been identified within IRDCs. 



Hot Molecular Cores (HMCs) are associated with the early stages of high-mass star 
formation. They correspond to the stage immediately after a cold dense starless core has 
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formed. During the later stages in their evolution, HMCs are often also associated with 
methanol masers and ultra-compact Hii regions. Numerou s examples of HMCs have been 



found throughout the Galaxy (e.g., iGarav fc Lizano 



1999 



Kurtz et al. 2000: Churchwell 



20021 ). including one within an IRDC fiRathborne et al.l 120071). HMCs are compact (<0.1 
pc), dense (~10^ — lO ^cm"^), and massive (~10^ Mrr^ jGaray fc Lizand Il999l : iKurtz et al. 



2000l : IChurchwelll 120021 ) . Due to the internal heating from the central high-mass protostar 
(T~100 K), they are typically very lurn inous (>10^ Lq) and show strong emission from com- 
plex molecules (e.g., Kurtz et al. 200oh. Accretion disks are als o inferred in th e later stage 



of the HMC phase direct l y ( ZhaneJ 



2000l : iBeuther et all l2002l : IZhang et al 



2005 



12001 



Cesaroni et al.l 120071 ) or indirectly (IKurtz et al. 



20051 ) by the presence of molecular outflows 



and maser emission. The existence of HMCs in IRDCs establishes a possible li nk between 
IRDCs and the early stages of high-mass star formation (IRathborne et al.l 120071 ). 



To understand the earliest stages in high-mass star formation one needs to identify and 
study the so-called 'high-mass starless cores,' the immediate precursors to high-mass proto- 
stars, HMCs, and ultra-compact Hii regions. High-mass star formation is rare and occurs 
rapidly; thus, the identification of high- mass starless cores is difficult. Because IRDCs have 
sizes, masses, and densities similar to cluster-forming molecular clumps but are considerably 
colder, we suggest that they are the precursors to clusters and their dense, cold cores the 
precursors to stars. Thus, the high-mass cores within IRDCs that contain no evidence for 
star formation are good candidates for the elusive 'high-mass starless cores.' 

Here we present interferometric sub-millimeter continuum and molecular line observa- 
tions obtained with the Submillimeter Array toward two cores within IRDCs; one showing 
evidence for active high-mass star formation (G034. 43-1-00. 24 MMl), the other showing none 
(G028. 53— 00.25 MMl). Both cores are at distances d > 3.7 kpc and h ave sizes R <0.8 pc , 
gas masses M >400 Mrr>, and bolometric luminosities L >10^ Lq (Table [T| ISimon et al.ll2006bl : 
Rathborne et al.ll2006l . 120081 ). We speculate that the 'active' core is in a later evolutionary 
stage than the more quiescent core. The new interferometric data reveal distinct properties 
for these cores and suggest that G034. 43-1-00. 24 MMl is a more evolved high-mass protostar 
in the HMC phase and that G028. 53— 00.25 MMl may be in a very early 'starless' core phase 
of high-mass star-formation. 
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2. Observations and Data Reduction 

2.1. Interferometric continuum and spectral line data 

Interferometric observations of G034.43+00 24 MMl and G028.53-00.25 MMl were 
carried out with the Submillimeter Arra}0 (SMA; Ho et al. 2004 ) on 2006 September 11 with 



seven antennas in the compact configuration. The projected baselines of the visibility data 
range from 15m to 70m. The double-sideband receivers cover a total of 4 GHz bandwidth. 
The digital correlator was configured to cover 329.1 to 331.1 GHz for the lower sideband and 
339.1 to 341.1 GHz for the upper sideband with a uniform channel spacing of 0.8125 MHz 
(~ 0.7km s^^) across the entire band. This setup resulted in simultaneous observations of 
the i^CO (3-2), C^^O (3-2), SO (73-67), and ^^SO (89-73) transitions. 

The 336 GHz zenith opacities, measured with the NRAO tipping radiometer located at 
the Caltech Submillimeter Observatory, were r ~ 0.4 (scaled from the 225 GHz measurement 
via T33GGHZ = 2.8 y<T225GHz)- The measured double-sideband system temperatures corrected 
to the top of the atmosphere were between 200 and 800 K. 

The primary beam of the SMA at these frequencies is ~ 37". The phase centers for the 
observations are hsted in Table [H We used the quasars 1751+096 and 1743—038 to cali- 
brate time dependent gains, and 3C454.3 to remove the gain variations across the passband. 
The fiux scale was referenced to Uranus. The visibility data were calibrated with the IDL 
superset MIR package developed for the Owens Valley Radio Observatory (OVRO) Inter- 
ferometer. The absolute fiux level is accurate to about 25%. After the calibration in MIR, 
the visibility data were exported to the MIRIAD format for further processing and imaging. 
The continuum is constructed from the line free channels in the visibility domain, and is 
further self-calibrated using the clean components of the image as input models. The gain 
solutions from the continuum self calibration is then applied to the spectral line data. The 
rms noise in the naturally weighted maps is ~ 12 and ~ 4 mJy beam~^ in the continuum (for 
G034.43+00.24 MMl and G028.53-00.25 MMl respectively). The noise in the continuum 
map for G034. 43+00. 24 MMl is higher because fewer line-free channels were available. For 
both cores the rms noise is 120 mJy beam^^ (~ 2.9 K) per 1.0 km s~^ channel in the line 
data. The final synthesized beam was ~ 2.3" x 2.1" which corresponds to a physical size 
of ~ 0.04 pc (~ 8,000 AU) and ~ 0.06 pc (~ 12,600 AU) for G034.43+00.24 MMl and 
G028.53-00.25 MMl respectively. 



^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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2.2. Single dish spectra 

To help with the interpretation of the SMA data, we have also included here lower- 
angular resolution spectra. These data were obtained as part of a large molecular line survey 
of cores within IRDCs conducted using the 30 m Institute de Radioastronomie Millimetrique 
(IRAM) telescope and the 15 m James Clerk Maxwell Telescope (JCMT). The molecular 
lines of CS (3-2), N2H+ (1-0), SiO (2-1), R^'^CO+ (1-0), and HCO+ (1-0) were obtained 
with IRAM in 2004 November, while the HCN (4-3), 1^00(3-2), and 0^^0(3-2) spectra 
were obtained with the JCMT in 2004 April. The IRAM spectra have a typical beam size of 
~ 26", an rms noise of 0.1 K channel"^, and a velocity resolution of 0.07km s~^. The JCMT 
spectra have a typical beam size of ~ 14", an rms noise of 0.1 K channel"^, and a velocity 
resolution of 1.1km s~^. 

In both cases, chopper wheel calibration, pointing, and focus checks were performed reg- 
ularly. All spectra were obtained in position-switched mode with a suitable nearby emission- 
free position used as a reference (typically ~5' from the core). The data were reduced using 
standard methods in the software packages GILDAS and SPECX. 



3. Results 

3.1. Low-angular resolution molecular line emission 

Figures [1] and [2] show the Spitzer 24 /um images of the two IRD Cs overlaid with the low - 



angular resolution (11") IRAM 30 m 1.2 mm continuum emission (iRathborne et al.ll2006l ). 
Note that while both IRDCs remain dark at 24 /xm, the core G034.43+00.24 MMl shows 
bright emission at 24 /im. 

Also included in these figures are the single dish low-angular resolution molecular 
line spectra obtained toward the cores. In both cases we see broad, saturated, and self- 
absorbed lines. In particular, toward G034. 43+00. 24 MMl we see broad line widths in 
the high-density tracers, CS (3-2) and HCN (4-3), shocked gas as traced by the bright 
SiO (2-1) emission, and high velocity line wings in HCO^ (1^0) and H^^CO^ (1^0)- To- 
ward G028.53-00.25 MMl we see saturated ^^CO (3-2) and N2H+ (1-0) line emission, faint 
emission from CS (3-2) and HCN (4-3) with broad line wings in the SiO (2-1), HCO+ (1-0), 
and H13C0+ (1-0) spectra. 

Table [2] gives a summary of Gaussian fits to each of these spectra: peak T^, central 
^LSR, and FWHM line width, AV. Because many of the lines show non-Gaussian profiles, 
the listed integrated intensity (/) was calculated by summing the emission over a suitable 
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velocity range [AV + 5 km s~^), rather than simply calculating the area under the Gaussian. 
For the N2H+ (1-0) spectra, however, we do not fit a single Gaussian profile, but instead 
fit the hyperfine components simultaneously. The derived opacity for the main hyperfine 
components are 0.10 for G034.43+00.24 MMl and 0.15 for G028.53-00.25 MMl. In these 
cases, the integrated intensity was calculated over the velocity range to include the emission 
from all hyperfine components. 

3.2. High-angular resolution continuum images 

Despite the similar sizes, masses, and bolometric luminosities (see Table [T]) derived for 
these cores from the low-angular resolution continuum images, the high-angular resolution 
sub-millimeter continuum images show significant differences (Fig. [3]). Toward G034. 43+00. 24 MMl 
the emission remains unresolved, while toward G028. 53— 00.25 MMl the emission is resolved 
into three individual conden^ationfl In both cases the emission is coincident with the peak 
in the lower-angular resolution millimeter continuum images (shown as contours on Figs. [T] 
and [2]). Table [3] lists the positions, deconvolved sizes, and peak fiuxes which were determined 
from Gaussian fits to the emission within the high-angular resolution images. 

3.3. High-angular resolution molecular line emission 

The high-angular resolution molecular line data toward these two cores are also very 
different. Toward G034. 43+00. 24 MMl the spectrum contains many molecular line emission 
features, but, in contrast, the spectrum toward G028. 53— 00.25 MMl contains no significant 
molecular line emission whatsoever. 

Figures m and O display the average spectra within a ~ 2" box centered on the peak in the 
high-angular resolution continuum image toward G034.43+00.24 MMl and G028. 53-00. 25 MMl 
respectively. Toward G034. 43+00. 24 MMl, the spectra show strong emission from the 
1=^00(3-2), Ci^O(3-2), 80(78-67), and ^^SOiSg-rg) transitions in addition to numerous 



^In this paper we use the term 'clump' to refer to IRDCs. These large (~ 5pc), massive (~ 10'^ M©), 
dense (~ lO'^cm"^) molecular structures are found within a Giant Molecular Cloud (GMC) and will likely 
give rise to a cluster of stars. Within these 'clumps', the many compact substructures are referred to as 
'cores' and have typical sizes of ~ 0.5 pc and masses of ~ 10^ Mq. At high-angular resolution, these 'cores' 
can often be further resolved into multiple substructures. We refer to the small continuum substructures 
identified by the SMA, with sizes ~ 0.05 pc, as 'condensations'. The 'condensations' presumably trace the 
material that will give rise to the individual stars. 
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lines from the complex molecules NH2CHO, CH2NH, CH3OCH3, CH3CN, and CH3CH2CN. 
Such complex molecules are typically produced and excited in the immediate surround- 
ings of a recently formed high-mass protostar as the grain ice mantles are evaporated (see 



Kurtz et al.ll2000l and references therein). Moreover, the presence of these species in such a 



compact region (< 0.03 pc) indicates that this is a HMC and, thus, an early stage in the 
formation of a high-mass protostar. 



4. Discussion 



4.1. Gas masses, Jeans masses, Virial masses, and Jeans lengths 

The optically thin sub-millimeter continuum emi ssion can be used to estimate the gas 
masses {Mq) of the condensations via the expression (IHildebrandl Il983l ) 



M, 



G 



where is the observed integrated source flux density, D is the distance, k,^ is the dust 
opacity per gram of dust, and BuiTa) is the Planck function at the dust temperature (Tp) . 



We estimate /«o.87mm by scaling the value of Ki.3mm (1-0 cm^ g~^: lOssenkopf fc Henninglll994l ) 
by v^, where /3, the dust emissivity index. We assume a gas-to-dust mass ratio of 100. In 
both cases we assume gray-body emission but with different dust temperat ures and emissiv- 



ity in dices. Fits to the IR-millimeter spectral energy distributions (SEDs; iRathborne et al. 



20081 ) for each core reveal dust temperatures of 34 K and 16 K and dust emissivity indices 
of 1.5 and 1.2 for G034.43+00.24 MMl and G028.53-00.25 MMl respectively (see Ta- 
ble [1]). We assume these values for f3 in the calculation of the condensation masses. For 
G034. 4 3-1-00. 24 MMl we assume a T/) of 100 K because this is a typical temperature for 
HMCs kurtz et al.lboooh . 



An assumption of the dust temperature for the condensations within G028. 53— 00.25 MMl 
is more difficult. Although fits to the core SED reveal a of 16 K, it is possible that the 
condensations may be small volumes of slightly warmer gas within this extended, cold enve- 
lope. Because strong emission at sub-millimeter wavelengths can trace either temperature or 
density enhancements ideally one needs high-angular resolution molecular line observations 
to directly measure the temperatures (e.g., NH3). The non-detection of any molecular lines 
in the high-angular resolution spectrum suggests that the temperature of the gas is quite 
cold. Thus, to calculate the masses for the condensations within G028. 53— 00.25 MMl we 
assume a temperature of 16 K. Because the temperature of the condensations may be higher 
than the temperature of the bulk of the gas within the core, the derived condensation masses 
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may, therefore, be an upper limit to the true mass. For example, if the condensations have 
a temperature of 30 K rather than 16 K, the derived masses will decrease by a factor of ~ 2. 

At high-angular resolution, the derived mass for G034. 43+00. 24 MMl is 29 Mq, while 
for G028. 53— 00.25 MMl the condensation masses are 9-21 (Table[3]). The lower- angular 
resolution single dish millimeter continuum images reveal these cores to have masses of 
430 and 1,130 respectiveljlf]. Thus, we find that only a small fraction (~ 4%) of the 
mass detected in the lower-angular resolution millimeter continuum data is recovered by the 
interferometer. Thus, it is likely that the emission traced by the interferometer corresponds 
to well-defined small, dense volumes of gas within a larger, more diffuse core. 

To determine the likelihood that the cores and condensations will fragment further, we 
calculate both their Jeans (Mj) and Virial (My) masses. These masses were calculated using 
the expressions 



Mj ^ 1.1 Mp, 





3/2 




10 K 




10 



p 



-19 g (.jj^ 3 



-1/2 



My = 1.145 X 10^ AV'^ DVn Mq 



where Tq is the dust temperature, p is the mass density (jZinnecker fc Yorkd 120071 ) . A V is 
the measured line width (km s~^), D is the distance (kpc), and f2 is the area on the sky 
(deg^). The mass density was estimated using the dust masses and by assuming the volume 
is equivalent to a sphere of radius R. 

The Jeans masses for the cores and condensations were calculated assuming the dust 
temperatures listed in Tables [1] and [3] and assuming that the gas and dust temperatures are 
identical. One must exercise caution, however, in applying a temperature of 34 K for the 
core G034. 43+00. 24 MMl as it is associated with a luminous 24 //m source which may be 
heating the gas to well above 34 K. For a gas temperature of 100 K, Mj will be higher by a 
factor of ~10. While the HMC's temperature may be closer to 100 K and correspond to a 
small volume in the core, the bulk of gas is likely to be colder. If the actual gas temperature 
is higher by a factor of two than the assumed value, then Mj will increase by a factor of ~4. 

In the case where thermal processes dominate, the ratio of the gas mass to Jeans mass 
(Mg/Mj) estimates the degree to which a gas structure is likely to gravitationally collapse 



■^These masses differ slightly from those quoted in iRathborne et al.l (|2006f ) because of their assumption 
of a single dust temperature (15 K) and emissivity index (2.0) for all their cores. 
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into thermally fragmented substructures. In contrast, if turbulence dominates, then the 
Jeans mass is no longer applicable and one must consider the Virial mass. In this case, the 
ratio of the gas mass to Virial mass (Mc/My) is a better indicator of further fragmentation. 
In either case, however, a ratio greater than unity indicates that the gas is unstable and will 
further fragment. 

Because the core line profiles in Figures [1] and [2] clearly show that the lines are much 
broader than what is expected from thermal broadening alone, it is likely that turbulence 
may significantly contribute to the pressure support within both cores. While we list the 
ratio Mg/Mj for the cores and condensations in Tables [1] and [3] it is unlikely that this ratio 
is physically significant. Instead, we use the ratio Mc/My to estimate the degree to which 
a gas structure is likely to gravitationally collapse into fragmented substructures. 

To calculate My for the cores we use the line widths measured from the low-angular 
resolution optically thin C^^O (3-2) emission (see Table[2]). For the condensations, we assume 
the same line width as measured toward the larger core. Because high-angular res olution NH.s 



obser vations toward a nearby IRDC show typical line widths of ~ 1.2 km s ^ (jWang et al. 



20061 ) on scales of the condensations, using the line width from the larger core may result in 
an overestimation of the calculated My for the condensations. 

We find that M^/My is 2.8 and 2.3 for the cores G034.43+00.24 MMl and G028.53-00.25 MMl 
respectively (Table [I]). For the condensations, we find that M^/My is < 1 (Table [3]). Since 
Mc/My is greater than unity for the cores, they are likely to fragment, but since M^/My is 
less than unity for the individual condensations, they are stable against further fragmentation 
and thus will probably give rise to individual stars. 



The Jeans length for the cores and condensations was calculated using the expression 



D 



Grupp 

where k is the Boltzmann constant. To is the dust temperature, G is the gravitational 
constant, nip is the mass of a proton, and p is the mass density. 

We find that the Jeans length, defined as the minimum length scale (radius) for grav- 
itational fragmentation to occur, is smaller (~0.04pc) for G034. 43-1-00. 24 MMl than that 
measured toward G028. 53— 00.25 MMl (~0.16pc; Table [1]). The separation of the con- 
densations within G028. 53— 00.25 MMl (~0.2pc) is comparable to the Jeans length of this 
core. Moreover, on the smallest scales of the individual condensations, the Jeans lengths are 
comparable to the measured radii of the condensations (Table [3]). These results support the 
idea that gravitational collapse is occurring within the IRDCs and that the fragmentation 
process has stopped at the smallest protostellar scales. 
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4.2. G034. 43+00. 24 MMl: a high-mass protostar with a rotating envelope? 



The core G034. 43+00. 24 MMl shows many indicators of high-mass star formation, in- 
cluding a bright 24 fim point sourc e, broad molecular line emission, shocked gas, water maser 



emiss i on, and molecular outflo ws (iGaray et al.ll2004j : IShepherd et al.ll2004l : iRathborne et al. 
20051 . bood : IWang et al.l bood ). G034.43+00.24 MMl also shows weak radio continuum 



emission consistent with a deeply embedded B2 protostar (IShepherd et al.ll2004l ). Fits to its 
IR-millimeter SED reveal it co ntains an object with a bolometric luminosity of ~ W^'^Lq, 
consistent with a BO protostar (IRathborne et al.ll2008l ). Because its high-angular resolution 
sub-millimeter continuum flux is dominated by a single, compact, massive emission feature, 
it is likely that the luminosity is dominated by a single high-mass protostar. 

Additional clues to the nature of this high-mass star-forming HMC can be obtained by 
examining the morphologies of its molecular line emission. For example. Figure [6] shows a 
distinct morphological difference between the CH3CH2CN and ^^CO emission. The 'hot 
core lines', in this case traced by the CH3CH2CN emission, are all unresolved at ~0.03pc 
resolution, while the ^^CO emission is easily resolved into a ~ 0.12 pc diameter extended 
region. Because the hot cores lines are excited in the immediate vicinity of a high-mass pro- 
tostar, it is likely, therefore, that the unresolved CH3CH2CN emission is tracing the heated, 
compact (<0.03 pc, < 6,500 AU) regions immediately surrounding the central protostar. On 
the other hand, the extended ^^CO emission is tracing the larger, less dense circumstellar 
region of the core. 

This extended ^^CO circumstellar emission also reveals a clear velocity gradient with 
red-shifted emission to the southeast and blue-shifted emission to the northwest of the central 
protostar (Fig. EJ. Such a velocity gradient might arise from either an outflow or a rotating 
envelope. Observational ev idence for accretion disks around high-mass protostars exists 
for only a few objects (see ICesaroni et al.l 120071 and references therein). The widespread 
existence of large (> 10,000 AU) disk-like structures in young HMCs would support the idea 
that high-mass stars form v ia monolithic collapse (e.g., iMcKee fc Taru 120031 ) rather than via 
competitive accretion (e.g. iBonnell fc Batdl2002l ). Because HMCs in IRDCs are probably 
very young, they are excellent candidates to search for high-mass accretion disks given their 
short lifetimes. 



In recent OVRO observations of CO (1-0) emission, IShepherd et al.l (120071 ) identify two 
outflows that may originate within the core G034. 43+00. 24 MMl. The dominant red-shifted 
lobe of the main outflow (their G34.4:A) appears to extend to the southwest, with little 
blue-shifted emission found to the northeast. Moreover, high-angular resolution SiO (1-0) 
imaging of this core with the VLA (Q. Zhang priv. comm.) confirms the presence and 
orientation of the CO outflow. This outflow appears to be orientated perpendicular to the 
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extended ^^CO (3-2) structure found here. Thus, the orientation of the outflow is consistent 
with this extended ^^CO (3-2) structure being a rotating envelope surrounding the central 
high-mass protostar. 

The position-velocity diagram of the circumstellar ^^CO (3-2) emission toward G034. 43+00. 24 MMl 
is consistent with Keplerian rotation around a Msin{i) ~ 10-50 M0 protostar (Fig. [7]). Thus, 
we see tentative evidence for a rotating circumstellar envelope surrounding this protostar. 
The derived mass for the central protostar is consistent with the luminosity and spectral 
type determined from the core's SED (10^'^ L©; BO). Higher angular and spectral resolution 
molecular line observations are required to better resolve the velocity field in order to confirm 
or refute the presence of an accretion disk. 

4.3. G028. 53— 00.25 MMl: an early phase of high-mass star-formation? 

In contrast to the bright, rich molecular line spectrum toward G034. 43+00. 24 MMl, 
the lack of strong molecular line emission toward G028. 53— 00.25 MMl suggests that it is 
either colder, more extended, or both. While ^^CO (3-2) and C^^O (3-2) emission is clearly 
seen toward this core in the lower-angular resolution data (Fig. [2l 14" beam), very little is 
seen in the high-angular resolution spectrum (Fig. \5\ upper limit of ~ 2.9 K). Moreover, in 
the 14" single dish spectra, G028.53-00.25 MMl shows only faint HCN (4-3) emission but 
does show self-absorbed CS (3-2), saturated or self-absorbed N2H+ (1-0), broad SiO (2-1) 
and H^^CO"^ (1^0)) ^ind an asymmetry in the HCO"'" (1-0) emission. The lack of emission 
from the high excitation HCN transition and the saturated N2H"'" profile suggests that the 
molecular gas is cold with a very high optical depth. Because SiO emission typically traces 
outflow activity caused by shocked gas, it is usually detected only in star-forming cores. 
Thus, the detection of broad line wings in the SiO, HCO"*" and H^^CO"*" spectra suggests 
that this core may already harbor some star formation activity. 

The absence of strong ^'^CO (3-2) or C^'^O (3-2) emission in the high-angular resolution 
spectrum may imply that the emission seen in the lower-angular resolution spectrum arises 
from a cold extended envelope which is resolved out by the interferometer. Using the limit 
of ~ 2.9 K for the brightness sensitivity of the SMA spectrum, we estimate from the single 
dish spectra that the sizes of the ^^CO (3-2) and C^^O (3-2) emitting regions must be > 16" 
(>0.4 pc). However, converting the derived mass to a ^^CO column density, we predict that 
the individual condensations have sufficient column densities (~ 10^^ cm~^) to be detected 
in ^^CO (3-2) emission in the current observations (peak T ~ 10 K and AV=5km s~^). It 
might be, therefore, that the ^^CO is depleted toward these condensation. It appears that 
G028. 53— 00.25 MMl does not contain any volume of heated gas and likely is a massive. 
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cold core. If a protostar does exist within G028. 53— 00.25 MMl then it has not yet heated 
its surrounding gas. Moreover, if it does contain a circumstellar disk, then it would be much 
bigger and colder compared to the disk surrounding G034. 43+00. 24 MMl. 

Given the core's mass is ~ 1000 M©, and assuming a star-formation efficiency of 30%, 
we estimate th at this core sh ould give rise to a cluster of 300 Mq. Using the expression rrimax 



1.2 M^i^^^^^ (lLarsonll2003l ). which relates the maximum stellar mass in a cluster 



JTimax) to 

the mass of the cluster {Mciuster)i "we estimate rrimax to be ~ 16 Mq. Thus, given its mass and 
a typical initial mass function (IMF), we would expect this core to give rise to ~ 4 high-mass 
(> 8 Mq) stars. Because it contains compact condensations, some evidence for star-formation 
activity but no obvious heating, we suggest that G028. 53— 00.25 MMl may be an example 
of a very early stage in high-mass star-formation. The measured bolometric luminosity 
(~1O^L0) may arise from the accretion of material onto the condensations. Given that the 
most massive condensation may further fragment it is possible that these cold condensations 
will continue to accrete material and give rise to a cluster of several high-mass protostars. 

It is interesting to compare our observations with recent models that suggest that the 
fragmentation of a core can be suppressed by the increase in its tempera ture due to accretion 
onto lower-mass protostars (jKrumholzll2006l : iKrumholz fc McKeell2008l ). In this scenario, a 
core's initial temperature and Jeans mass predicts it will fragment into many low-mass con- 
densations. However, as these condensations form and begin to accrete, the temperature 
increases which, in turn, increases the overall Jeans mass of the core. An increase in tem- 
perature to ~ 100 K is thought to be sufficient to halt the fragmentat i on pro cess. To form 
high-mass cores and avoid further fragmentation, iKrumholz fc McKed (120081 ) find that the 



cloud's critical column density, E, should be > Igcm ^ and its luminosity-to-mass ratio, 
L/M, > 10 Lq/Mq. For the core G028.53-00.25 MMl we find that S is ~ 0.4 gcm'^ and 
L/M is ~ 0.3 Lq/Mq (Tabled]). These rough estimates are smaller than the predicted thresh- 
olds. Thus, fragmentation and the formation of lower-mass protostars is still possible within 
G028.53-00.25 MMl. This is in contrast to G034.43+00.24 MMl which has a measured S 
and L/M much greater than these thresholds implying that it will give rise to a high-mass 
star. In this case, however, its high Mg- to Mj ratio (for Id = 100 K) suggests that the 
increase in gas temperature alone may not be enough to halt thermal fragmentation. 

The ongoing fragmentation within G028. 53— 00.25 MMl is consistent with our detection 
of multiple condensations toward its center. If this model is correct, then it appears that 
only the condensation A has sufficient column density to form a high-mass star, the other 
condensations will most likely only form low-mass stars (Table |3]). The derived temperature 
for the dust within this core is well below 100 K, suggesting that material has not been heated 
sufficiently to halt the fragmentation. Due to its high mass, this represents an excellent core 
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in which to test these theoretical predictions for high-mass star-formation. 



5. Conclusions 

Using the SMA we have obtained high-angular resolution sub-millimeter continuum im- 
ages and molecular line spectra toward two cores within IRDCs. In the low- angular resolution 
data these cores have sizes R <0.8pc, g as masses M >4OOM0, and bolometric luminosities 



L>102Lq flRathborne et al.ll200d 120081 ). Despite their similar sizes, masses, and bolometric 
luminosities the two cores show remarkable differences in the high-angular resolution data. 
The sub-millimeter continuum images reveal that G034. 43+00. 24 MMl remains unresolved 
and has a size of ~ 0.03 pc and mass of 29 Mq. In contrast, G028. 53— 00.25 MMl contains 
at least three compact (~ 0.06 pc), massive (9-21 M0) condensations. 

Moreover, the molecular line spectra reveal that G034. 43+00. 24 MMl emits strong lines 
from many complex molecules, but the core G028. 53— 00.25 MMl does not. This suggests 
that G034.43+00.24 MMl is hot, compact, and dense, while G028. 53-00.25 MMl is colder 
and possibly more extended. We speculate, therefore, that G034. 43+00. 24 MMl is a HMC 
in an early stage in the formation of a high-mass protostar, while G028. 53— 00.25 MMl 
may represent an even earlier phase in the high-mass star-formation process. Because 
G028. 53— 00.25 MMl shows some evidence for active star formation in the low-angular reso- 
lution data (e.g., shocked gas, high velocity line profiles, and a high bolometric luminosity), 
we cannot rule out the possibility that the condensations within this core may have already 
begun to accrete. However, the low temperatures and absence of bright 24 /jm point sources 
or strong molecular line emission from high excitation transitions suggests that there is little 
localized heating. Because of its size and mass, G028. 53— 00.25 MMl is likely to give rise 
to at least one high-mass star. We speculate that the detected condensations may continue 
to accrete material from their surroundings to eventually form high-mass stars. Since the 
condensations account for only ~4% of the total mass, there is an ample amount of supply 
of material to accrete. 

The properties of these cores support the idea that high-mass star and cluster forma- 
tion likely occurs in IRDCs. High-angular sub-millimeter continuum images and molecular 
line spectra toward more cores within IRDCs are needed to unambiguously establish the 
evolutionary path from a high-mass starless core to a high-mass protostar. 
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Table 1: Core properties 



Core 


Coordinates 




D 


R 


Td 


/3 


L 


Mg 


Mg 
Mj 


Mg 
My 


Aj 


S 


L 
M 




RA 


Dec 


























(J2000) 


(J2000) 


(km s-i) 


(kpc) 


(pc) 


(K) 




(L©) 


(Mq) 






(pc) 






G034.43+00.24 MMl 


18:53:18.0 


01:25:24 


57 


3.7 


0.10 


34 


1.5 


19,600 


430 


550 


2.8 


0.04 


3.2 


46 


G028.53-00.25 MMl 


18:44:18.0 


-03:59:34 


87 


5.7 


0.42 


16 


1.2 


370 


1130 


800 


2.3 


0.16 


0.4 


0.3 


References 






1 


1 


2 


3 


3 


3 


3 


4 


4 


4 


4 


4 



References, 
work. 



(1) ISimon et all |2006bh : (2) iRathborne et all |2006l ): (3) iRathborne et al.l |2008l ): (4) this 



Table 2: Summary of the low- angular resolution molecular line spectra 



Molecular line 




G034.43+00.24 MMl 




G028.53 


-00.25 MMl 








AV 


ja 


'-A 




AV 






(K) 


(km s"-') 


(km s-i) 


(Kkm s-l) 


(K) 


(km s-i) 


(km s-i) 


(K km s-i) 


ClSO(3-2) 


2.9 


57.8 


3.7 


11.1 


1.2 


87.0 


3.2 


4.8 


13C0 (3-2) 


8.2 


57.6 


5.1 


44.8 


2.4 


86.2 


5.8 


14.8 


HCN (4-3) 


3.4 


57.7 


12.5 


45.8 


1.2 


86.0 


7.6 


1.1 


CS (3-2) 


4.1 


57.9 


6.8 


31.7 


0.5 


86.1 


7.4 


4.2 


N2H+ (1-0) 


l.O'' 


57.8 


2.6 


57.1 


0.2*= 


86.8 


3.6 


16.4 


SiO (2-1) 


0.8 


57.9 


10.6 


9.3 


0.3 


86.4 


9.0 


3.2 


H"cO+ (1-0) 


0.8 


57.9 


10.5 


9.0 


0.3 


86.4 


9.1 


3.3 


HCO+ (1-0) 


2.5 


60.2 


3.7 


13.4 


1.1 


87.5 


9.2 


11.2 



"The integrated intensity quoted here is not the area under the Gaussian, but instead is calculated over a 
velocity range to include all the emission. 
''This is tT^ for the main hypcrfine components. 



Table 3: Condensation properties 



Condensation 


Coordinates 
RA Dec 


Peak 
flux 


Deconvoled 
diameter 


Td 


/3 


Mg 


Mg 
Mj 


Mg 
Mv 


Aj 


s 




(J2000) 


(J2000) 


(mJy) 


(pc) 


(K) 




(Mq) 






(pc) 




G034.43+00.24 MMl 


18:53:18.03 


1:25:25.58 


2760 


0.03 


100 


1.5 


29 


30 


1.2 


0.02 


7.7 


G028. 53-00.25 MMl A 


18:44:18.08 


-3:59:34.33 


73 


0.05 


16 


1.2 


21 


160 


0.8 


0.01 


2.7 


G028. 53-00.25 MMl B 


18:44:17.95 


-3:59:31.27 


33 


0.10 


16 


1.2 


9 


10 


0.2 


0.07 


0.2 


G028.53-00.25 MMl C 


18:44:17.93 


-3:59:35.42 


31 


0.06 


16 


1.2 


9 


30 


0.3 


0.03 


0.8 
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Fig. 1— G034.43+00.24. Left: Syitzer 24 urn im age with contours of the IRAM 30 m 



1.2 mm continuum emission (IRathborne et al. 



20061 : contour levels are 60, 90, 120, 240, 480, 



1200, and 2200 mJy beam~^; color scale is logarithmic from 30 MJy sr"^ [black] to 200 MJy 
sr~^ [white]). Note that while the IRDC remains dark at 2Afim the core MMl is associated 
with bright 24 fim emission. Right IRAM and JCMT molecular line emission toward the 
core MMl. The solid vertical line marks the central velocity of the core, as traced by the 
optically thin C^^O (3-2) emission. Note the broad line widths in the high-density tracers, 
CS (3-2) and HON (4-3), the shocked gas as traced by the bright SiO (2-1) emission and 
the high velocity line wings in HCO"*" (1-0) and H^^CO"'" (1-0). Table [2] summarizes the 
results of fits to these spectra. 
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Fig. 2.— G028.53-00.25. Left: Syitzer 24 urn im age with contours of the IRAM 30 m 
1.2 mm continuum emission ( IRathborne et al.ll2006l : contour levels are 60, 90, 120, 180, 240, 
and 320 mJy beam~^; color scale is logarithmic with 50 MJy sr^^ [black] to 100 MJy sr~^ 
[white]). Note that the IRDC and, in particular, the core MMl remain dark at 24 /im. Right 
IRAM and JCMT molecular line emission toward the core MMl. The solid vertical line 
marks the central velocity of the core, as traced by the optically thin C^^O (3-2) emission. 
Note the saturated ^^CO (3-2) and N2H+ (1-0) line emission and the faint emission from 
CS (3-2) and HON (4-3). The presence of broad line wings in the SiO (2-1), HCO+ (1-0), 
and H^'^CO^ (1^0); spectra suggests that some star-formation activity is occurring within 
this core. Table [2] summarizes the results of fits to these spectra. 



- 20 - 



o 



-10' 



G034. 43+00. 24 MMl 



- O Beam 



0.04 



I 



c 



-3°59'20" ^ 



1000 



iB^sa^ig.o^ ia.5= ib.o=' 17.5= 
Right Ascension (J2000) 




140 



18"44™19.0= 18.5= 16.0= 17.5= 
Riglit Ascension (J2000) 



Fig. 3. — SMA continuum images toward the two cores (the images cover ~ 30" around 
the core's center). The contour levels for G034.43+00.24 MMl (left) are 0.3 to 3.0 in 
steps of 0.5 Jy beam^^ and for G028.53-00.25 MMl (right) are 0.015 to 0.07 in steps of 
0.01 Jy beam~^. In both cases the lowest two contour levels are in black for clarity. Note 
that G034. 43+00. 24 MMl contains a single, compact (~ 0.03 pc), massive (24 M©) emission 
feature, while G028. 53— 00.25 MMl appears resolved into three individual condensations. 
These condensations are marked with a cross and have sizes of ~ 0.06 pc and masses of 
8-19 Mq. Table [3] list their properties. 



- 21 - 



15 F"^ 




I I I I , I I I I I I I I I I I 

329.5 330 330.5 331 

Frequency (GHz) 



Fig. 4.— SMA spectra toward G034.43+00.24 MMl. Both the lower and upper sidebands 
are shown (upper and lower panels respectively). The spectra are averaged over a ~ 2" 
region centered on the peak in the high-angular resolution continuum image. Note the strong 
emission from the many complex molecules. Thus, we suggest that G034. 43+00. 24 MMl is 
a HMC associated with the very earliest stage in the formation of a high-mass protostar. 
The molecular lines marked on this figure were based on the output from the Lovas on-line 
catalog of suggested rest frequencies. 
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Fig. 5. — SMA spectra toward G028. 53— 00.25 MMl. Both the lower and upper sidebands are 
shown (upper and lower panels respectively). The spectra are averaged over a ~ 2" region 
centered on the peak in the high-angular resolution continuum image. While ^^CO (3-2) 
and C^^O (3-2) emission is clearly seen toward this core in the lower-angular resolution data 
(Fig. [2j 14" beam), very little is seen in this spectrum. The absence of strong ^^CO (3-2) or 
C^^O (3-2) emission in the high-angular resolution spectrum may imply that the emission 
seen in the lower-angular resolution spectrum arises from a cold extended envelope which is 
resolved out by the interferometer. 




1B''53°1B.40=1B,20= IB. 00' 17,80= 17,60' 
Right Ascension (J2000) 

Fig. 6. — SMA CH3CH2CN integrated intensity image in gray scale tracing the unresolved 
HMC within G034.43+00.24 MMl. Contours show the red-shifted {Vlsr= 61km s"^) and 
blue-shifted (yLSR= 55km s^^) ^^CO (3-2) emission surrounding the unresolved, bright 
HMC (emission is from a single channel, with the same contour levels for the red and blue 
components; 1.7 to 7.5 in steps of 0.5 Jy beam~^). This emission appears to be tracing an 
extended rotating envelope. 
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Fig. 7.— SMA i3C0(3-2) position-velocity diagram for G034.43+00.24 MMl. The slice 
is taken perpendicular to the axis of the rotating envelope and shows the two components. 
Overlaid on this image are models for Keplerian rotation with central masses of 10, 20, 50, 
and 100 Mq. The emission from G034. 43+00. 24 MMl is consistent with Keplerian rotation 
around a ~ 10-50 source. 



